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ABSTRACT 

Cosmic rays are appealing as a source of ionization in starburst galaxies because of the great columns they 
can penetrate, but in the densest regions of starbursts, they may be stopped by pion production and ionization 
energy losses. I argue that gamma rays are the source of ionization in the deepest molecular clouds of dense 
starbursts, creating Gamma-Ray Dominated Regions (GRDRs). Gamma rays are not deflected by magnetic 
fields, have a luminosity up to ~ 1 /3 that of the injected cosmic rays, and can easily penetrate column depths 
of ~ 100 g cm" 2 before being attenuated by 7Z pair production. The ionization rates of GRDRs, < 10~ 16 sec -1 , 
are much smaller than in cosmic ray dominated regions, but in the most extreme starbursts, they may still 
reach values comparable to those in Milky Way molecular clouds. The gas temperatures in GRDRs could be 
likewise low, < 10 K if there is no additional heating from dust or turbulence, while at high densities, the kinetic 
temperature will approach the dust temperature. The ratio of ambipolar diffusion time to free-fall time inside 
GRDRs in dense starbursts is expected to be similar to those in Milky Way cores, suggesting star-formation 
can proceed normally in them. The high columns of GRDRs may be opaque even to millimeter wavelengths, 
complicating direct studies of them, but I argue that they could appear as molecular line shadows in nearby 
starbursts with ALMA. Since GRDRs are cold, their Jeans masses are not large, so that star-formation in 
GRDRs may have a normal or even bottom-heavy initial mass function. 

Subject headings: galaxies: starburst - galaxiesTSM - gamma rays: galaxies - cosmic rays - galaxies: star 
formation 



1. INTRODUCTION 

With large star-formation rates often concentrated into 
small volumes, starburst galaxies are filled with the many 
kinds of high energy radiation that accompany star-formation, 
from ultraviolet photons to cosmic rays (CRs). This radia- 
tion can interact with the interstellar medium (ISM) of the 
abundant gas to heat or ionize it. The ionization of gas in 
star-forming environments in particular can drive chemistry 
in molecular clouds and sets the dynamics of magnetic fields 
that may in turn regulate star formation. 

However, the ionization rate from this radiation is limited 
by two key factors: the luminosity of the radiation, and the 
regions within the starburst it can penetrate into. Energy con- 
servation ultimately limits the rate atoms can be ionized by 
radiation; more radiation generated means more gas can be 
ionized, subject to the efficiency with which that radiation is 
converted into ionization. A powerful source of ionizing lumi- 
nosity are extreme ultraviolet (EUV; hv > 13.6 eV) or Lyman 
continuum photons from young O stars. The O stars create H 
II regions around themselves, where the gas is highly ionized. 
Indeed, the ionizing photon rate as measured by free-free ra- 
dio luminosity is correlated w ith the star-formation rate (e.g., 
ICondonl \\99% iMurphviHoTlh . The amount of ionized ma- 
terial in H II regions is large because of the high efficiency 
that EUV light ionizes surrounding gas: starburst galaxies are 
highly optically thick to EUV light (e.g., Leithere r et al J 1 995t 
iHurwitz et aDl!997t iHeckman et alJl2001l) . Almost all of the 
power in EUV radiation goes into ionizing H II regions or 
possibly into dust absorption. 

However, the high optical depth of neutral gas to EUV light 
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also poses a problem for the ionization of most of the ISM. 
Most of the EUV light is downgraded before it goes very far, 
so that H II regions are a small fraction of the volume of star- 
burst galaxies, and something else must be responsible for 
ionizing the dense molecular gas in starburst galaxies. Far 
ultraviolet (FUV; 6 eV < hv < 13.6 eV) light is not absorbed 
by neutral hydrogen, but it can dissociate molecular hydro- 
gen and ionize carbon: thus, regions where FUV light pen- 
etrates into are p hotodissociation regions (PD Rs; see the ex- 
tensive review by Hollenbach & Tielens 1999). PDRs include 
much of the neutral gas in the Milky Way, but FUV light is 
quickly absorbed by dust (limiting both the power going into 
ionization and the penetration depth), and PDRs extend only 
to columns of Nh ~ 2 x 10 22 cm" 2 (0.03 g cm" 2 ), smaller than 
typical columns through st arbursts (0 . 1 - 1 g cm" 2 ; iKennicutl 
19981 iHopkins et al.ll2010l) . 

The interiors of molecular clouds are ionized by more pen- 
etrative but less luminous sources of radiation. X-rays can 
traverse higher columns and can ioniz e neutral gas, creat- 
ing X-ray Dom inated Regions (XDRs; Maloney_etal][l996; 
Meiierink & Spaans 2005). In active galactic nuclei with high 
X-ray luminosities, X-rays can provide most of the ioniz- 
ing lu minosity, although t his is concentrated towards the nu- 
cleus (Papadopoulos 2010). Even in pure starburst galaxies, 
hard X-rays are emitted by high mass X-ray binaries, with 
a typical luminosity of Lx ~ lO -4 !*, where L + is the bolo- 
metric (largely in the infrared) luminosity fro m star-formation 
(e.g.. iRanalli et al.ll2003t iPersic et al.ll2004l) . However, most 
galaxies are optically thin to hard X-rays, meaning most 
of that luminosity simply leaves the starburst before ioniz- 
ing neutral gas. Furthermore, even hard X-rays are stopped 
once the columns become Compton thick {Nh ~ 10 24 cm" 2 ; 
2.5 g cm" 2 ), which can o ccur in extreme starbursts like Arp 
220 dLehmer et al.l [2010) and in molecular clouds. Even 
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higher column densities can be attained in the parse c-scale gas 
around active galactic nuclei ( Hop kins et al.ll20l2l and refer- 
ences therein). 

An attractive candidate for the ionization of starbursts is 
cosmic rays 0, which ionize the interiors of Milky Way 
molecular clouds. The CRs regulate the temperature and 
ionization fraction in molecular clouds and cores, and be- 
cause they pervade Milky Way molecular clouds, they 
set uniform conditions for star-formation throughout the 
Galaxy. CRs have many advantages as ionization agents 
in starbursts. (1) They are produced by star-formation 
at similar luminosities as X-rays, with Lcr ~ 3 x lO" 4 /^ 
(iThompson. Ouataert. & Waxmanl 120071) . Since they trace 
star-formation, they naturally should be present in star- 
forming regions, unlike X-rays from AGNs, which are intense 
only near the nucleus itself. (2) They can traverse columns 
of up to ~ 80 g cm" 2 at GeV energies, stopped mainly by 
inelastic collisions with ISM atoms that produce pions. At 
lower kinetic energies (less than a few hundred MeV), ion- 
ization losses also stop them, becoming rapidly more effec- 
tive at smaller kinetic energies. They can therefore penetrate 
into regions of high dust extinctions, a big advantage over UV 
photons in the dusty environments of starbursts. (3) Galaxies 
are extreme scattering atmospheres for CRs, deflecting them 
with magnetic inhomogeneities. This allows the CR energy 
densi ty to buil d up to very high l evels (e.g.. ISocrates et all 
2008). Thus, Papadopoulos (2010) proposed that molecular 
gas in starburst g alaxies are Cosmic Ra y Dominated Regions 
(CRDRs; see also lSuchkov et alJI 19931) . 

But the high scattering optical depth to CRs is a double- 
edged sword. On the one hand, it means that each CR has 
many chances to ionize gas atoms as it scatters across a star- 
burst, but on the other it means that a CR has many chances to 
be destroyed through pionic and ionization losses. The scat- 
tering atmosphere means that CRs see a much higher col- 
umn that can stop them than photons do. Indeed, in the 
Milky Way, there is evidence that low energy (< 100 MeV) 
CR components contribute to high ionization rates in the 
diffuse ISM but are stopped by their own ionization losses 
dlndriolo et al.l [2009b . In dense starburst galaxies, protons 
may even lose all of their energy to pionic (and ionization) 
losses, in the "proton cal orimeter" l imit - with little of the 
CR luminosity escaping dPohlll 1 994t lLoeb & Waxmanll2006t 
Thompson, Oua taert. & Waxmanl 120071: lLackietal.1 120 10i ). 
So while starbursts with their dense gas may efficiently con- 
vert CR energy into ionization, CRs may be unable to reach 
heavily shielded molecular regions in starbursts. The GeV 
and TeV gamma-ray observations of M82 and NGC 253 
(lAcero et al]l200% lAcciari et al.ll2009T: lAbdo et al.ll2010l) in- 
dicate that roughly f w ~ 20 - 40% of their CR power is lost 
to pionic emission ( Lac k feTaTI 1201 ll) . much higher than in 
the Milky Way w here the ratio is only a few percent (e.g., 
Stron g~etaTI 120101) . implying that pionic losse s may be im- 
portant for CRs in starbursts (Lacki et al. 20Tob. 

An even more worrying problem wit h CR ionization in 
starbursts are the pres ence of winds (e.g., Cheva lier & Cleggl 
U985tlHe ckman 2003 1): the se carry CRs out of their starbursts 
(c.f., ISuchkov et al.l 119931) . effectively wasting the power 
they could provide in ionizing the ISM; but unlike diffu- 

3 In genera l, when I refer to CRs in this paper , I mean CR protons (and 
other nuclei) ISuchkov et al. 1993; Papadopoulos 2010), which make up the 
bulk of the CR power in the Milky Way. However, low energy CR electrons 
and positrons can also contribute to the CR ionization rate. 
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FIG. 1. — Schematic diagram showing the kinds of ionizing radiation in 
star-forming galaxies, and the regions they can ionize. UV light has a great 
luminosity, but cannot penetrate far into the starburst; gamma rays are low 
luminosity but can penetrate deep into molecular clouds. Cosmic rays and X- 
rays from X-ray binaries are intermediate in both luminosity and penetration 
depth. The figure is not to scale, and the units are arbitrary. 



sion, wind s do not let them p enetrate deeper into molecu- 
lar clouds (ICrocker et al.l 1201 1[). The hard gamma-ray spec- 
tra of M82 and N GC 253 dAce ro et al l 120091 lAcciari et al.l 
120091; lAbdo et all 120101) . combined with the evidence that 
fx < 1, suggest that winds are important for CR trans- 
port in them, though denser starbur sts like Arp 220 may be 
true proton calorime ters. However, Suchk ov et al.lll993l and 
Papadopoubi|20i3 did assume that winds limited the life- 
times of CRs in starbursts when calculating the CR energy 
densities. 

To summarize the problem for ionizing radiation, radia- 
tion is optimally effective in uniformly ionizing regions when 
r = 1 : in media with smaller r, the radiation escapes with- 
out efficiently ionizing the gas, but in media with larger r, 
the radiation ionizes only the skin of the region and leaves a 
shielded core. We therefore expect a hierarchy of ionization 
regions in starburst galaxies, with each step having lower lu- 
minosity but higher penetration depth. The H II regions and 
PDRs from UV light receive the most ionizing power but only 
constitute a small volume of starbursts; XDRs and CRDRs re- 
ceive much less power but pervade denser phases of starbursts 
(Figure [TJ. Is there anything beyond CRDRs, which could 
guarantee the ionization of the molecular cores of starbursts? 

I argue here that this ionization could come from gamma 
rays. High energy gamma rays are produced in starbursts 
by the pionic losses of protons. Additional leptonic gamma 
rays come from electrons and positrons (e ), some of which 
may be secondary made by pion production. Unlike CRs, 
gamma rays are not deflected by magnetic fields, and they are 
stopped only by Bethe-Heitler 7Z pair production at columns 
of ~ 10 26 cm" 2 (~ 160 g cm" 2 ). In 7Z pair production, the 
gamma ray produces an electron and a positron in the elec- 
tric field of an atomic nucleus (Z+7 — >Z + e~ + e + ) (e.g., 
iBerest etskii et a l.l 19791) . The result in dense columns is a cas- 
cade: the either ionize the ISM directly if they are low 
enough energy, or effectively downgrade the energy of the 
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initial gamma ray when they cool by radiating. Thus it is 
possible to use the proton calorimetry of a starburst against it: 
the gamma rays can act as a "secon d wind" for CR ionization 
(c.f., [Omebayashi & Nakano 1981). The result is that even if 
CRs are stopped in the outer reaches of a molecular cloud, the 
inner regions can still be ionized as a Gamma-Ray Dominated 
Region (GRDR; FigureQ]). Of course, in most of the starburst, 
gamma rays escape most of the time, with most of the energy 
"wasted" that way, but the key point is that gamma rays pro- 
vide a low level of ionization that extends extremely deep into 
the molecular clouds of a starburst. 

This paper will consider basic properties of GRDRs in 
extreme starburst environments. Using the energetics of 
gamma-ray production in starbursts, I estimates the charac- 
teristics ofGRDRs in section |2j including the ionization rate 
(section l2~TT) , ionization fraction and am bipo lar diffusion time 
(section l2~3l l. and temperature (section |2~4V I then estimate 
when the transition from CR to gamma-ray ionization occurs 
and demonstrate how the extreme scattering atmospheres of 
starburst can turn even clouds with relatively small columns 
into GRDRs (section [3]). To more accurately calculate how 
the pair e ± ionize the ISM, I construct a simple ID model 
of a 7Z cascade in a starburst in section @] I consider how 
GRDRs may appear observationally, including their indirect 
effects on stellar initial mass functions in section [5] I finally 
discuss some additional speculations in the Conclusions, sec- 
tion |6] 

2. BASIC CHARACTERISTICS OF GRDRs 
2.1. Estimates of GRDR Ionization Rates 

Suppose a starburst is a disk with radius R and a midplane- 
to-edge scale height h. It has a gamma ray luminosity of L 7 , 
giving a gamma-ray energy flux of F 7 = Lj/(2ttR 2 ). The vol- 
umetric ionization power that gamma rays deposit in a gas 
cloud that is optically thin to gamma rays (r 7 z < 1) is 



■1 

ion 



(i) 



where I is the gas column length, Hh is the number density of 
hydrogen atomfl, Nh = is the column density, <r 7 z is the 
cross section for 7Z pair production, and f^ on is the fraction 
of energy in pair e ± that goes into ionization losses. The f? on 
quantity depends on the physical conditions within the cloud 
and the spectrum of the gamma rays; here I will estimate it as 

/I = 0.1. 

When a typical high energy proton ionizes an atom, the 
ejecte d electron will have a typica l kinetic energy of ~ 35 eV 
(e.g., ICravens & Dalgarriol 119781) : the CR proton therefore 
typically loses (35 + 1 3 .6) = 49 e V per ionization event. On the 
other hand, the initial ionized electron i n turn ionizes, on av- 
erage , 4> = 0.7 additional electrons (e.g., Cra vens & Da lgamo 
1978). Therefore, if ionizations by high energy e ± are sim- 
ilar to those by high energy protons, each e^ requires E\ on = 
49 eV/(l + </>) = 29 eV per ionization event, whether primary 
or secondary. Here I scale Ei on to 30 eV for simplicity. 

The ionization rate from gamma rays is then simply £ 7 » 
rlJ(n H E i0B ): 

C 7 «^5^. (2) 



4 Throughout this paper, I use the number density of hydrogen atoms, hh, 
not the number density of hydrogen molecules niHj). The conversion be- 
tween the two is riH = 2n(/?2). Likewise, the column density Nh is for hydro- 
gen atoms, not hydrogen molecules. 



Few starbursts are gamma-ray detected, but we can estimate 
the gamma-ray flux using the star-formation rate, the Schmidt 
law, and the time scale for pionic energy losses. The gamma- 
ray luminosity is 



(3) 



where f n is the fraction of the CR power that goes into pionic 
losses: /„. w t^ & /t v = [1 + ^Awm] -1 at high energies. The 
pionic loss time is roughly 



i 50 Myr 



cm 



(4) 



where (ng ) is the mean density of gas in the entire starburst 
dMannheim & S chlickeiser 1994); this is equivalent to saying 
pionic losses stop CR protons over a column = t^nncmH = 
79 g cm" 2 . In terms of the mean column density of gas in the 
starburst (S g ) = 2m.nhnu, we have t n ~ 100 Myr mnh/ '(£#)■ 
The wind escape time is f w i nc j = /z/v w i n d- Combining these 
timescales, I have: 



1+0.16 



"wind 



300 km s- 



g cm _: 
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The luminosity of injected CRs scales with the su- 
pernova rate (or more generally, the massive star for- 
mation rate). A CR luminosity of 10 50 erg per su- 
pernova is expected. For a Salpeter IMF, we expect 
a supernova rate of r SN = 0.0 64 y r^SFR/Mg yr" 1 ) 
(Thompson. Ouatae rt. & Waxmanl 120071) . In turn, I 
can relate the star-formation rate to the gas den- 
sity using the Schmidt law, Esfr = SFR/(ttR 2 ) = 
2.5 x IO^Mq yr" 1 kp C - 2 ({E g )/M pc" 2 ) (lKennicuttlll99^) . 
Putting these relations all together, I find the gamma-ray 
flux: 



F 1 s» 0.0036 erg sec cm f„ 



F 1 « 0.0252 erg sec cm / w 



( R 



SFR 

M yr" 1 J V250pc 



gem 



(6) 
(7) 



and finally, 

0.157 ergsec-'em" 2 (g^r) 24 ((E,) < 0.16 g cm" 2 ) 

0.0252ergsec- I cm- 2 (g^ i y ' 4 ((S g ) > 0.16 g cm" 2 ) 

(8) 

Evaluating the cross section for 7Z pair production 
dBerestetskii et al.ll 19791) . 
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at £ 7 = 1 GeV to be 1 .02 x 10" 26 cm 2 , 1 find 
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FIG. 2. — Plot of ionization rate from gamma rays (solid) and CRs (dotted: 
momentum power law injection; dashed: total energy power law injection), 
using the Schmidt Law. The GRDR ionization rate is much smaller than the 
CRDR ionization rate at all surface densities, but the ratio of £ 7 / fcR grows 
with (E g ). Gamma-ray ionization is more important than 40 K radioactivity in 
virtually all starbursts, and approaches the £ value in Milky Way molecular 
clouds at the highest gas surface densities. Individual galaxies may have 
higher or lower f depending on their deviation from the Schmidt Law. I have 
assumed v w ind = 300 km s and =0.1. 

The results are plotted in Figure [2] These ionization rates are 
small, though they do rise to ~ 10~ 17 sec -1 for the densest 
ULIRGs with w 10 g cm" 2 . In principle, if /; on is 1, 
these ionization rates can be ten times higher. In any case, 
though, they remain much higher than the ionization rate from 
long-lived radioactive elements such as 40 K in the Milky Way, 
Cr « 10" 22 sec -1 (lUmebavashi & Nakanolll98lL 120091) even 
for mean column densities of w 0.01 g cm -2 . 

In starbursts, there may be additional radioactive ionization 
from relatively short-live d isotopes like 26 A l because of the 
high supernova rate (c.f., Diehl et al. 2006); for example, in 
the young Solar System, 26 Al alone sustained an ionization 
rate of ~ 10~ 18 sec -1 (lUmebavashi & Nakano 2009; see also 
Step inskil[T992h . In a companion paper (Lacki 20120, I find 
that the mean 26 Al ionization rate in starbursts quite high be- 
cause of their rapid specific star-formation rates, of the order 
10~ 18 -10~ 17 sec" 1 . If true, gamma-ray ionization only dom- 
inates in the most extreme starbursts with ss 10 g cm" 2 . 
However, whether these high ionization rates are actually at- 
tained in the molecular gas of starbursts depends on whether 
the 26 Al is actually well-mixed with the gas; as with CRs, the 
propagation of 26 Al from its sources matters. This is not an 
issue with the gamma rays. 

Some starbursts are in fact gamma-ray detected, allowing 
us to directly calculate F~. The brightest starburst in the 
gamma-ray sky in terms of observed flux is M82, with a GeV- 
TeV g amma-ray luminosity of 2 x 10 40 e rg s" 1 (|Abdo et all 
2010) and a starburst radius of 250 pc (e.g. JGoetz et al.lll990b 
Williams & Bowerl HoiO). Therefore, using eqn.[2] I find 

C7 = uxir.W.(^)(J^f, (,„ 



much smaller than even in Milky Way molecular clouds, but 
still much larger than the ionization rate from long-lived ra- 
dioisotopes like 4() K. 

Arp 220, which if star-formation powered has a star- 
formation rate of 200-300 M^ vr" 1 dTorre si 120041 and ref- 
erences therein: lAnan tharamai ah et al . 2000), is a more likely 
home for GRDRs. Most of its radio emission comes from two 
50 - 1 00 pc radius nuclei (iNorrisll 19881 iDownes & Solomonl 
1998), although there are varying estimates of the literature 
of how much of the infrared luminosity is actually from the 
nuclei (IDownes & Solomonl 119981; IDownes & Eckaril 120071; 
ISakamoto et al.ll2008h . If I assume that each nucleus has a 
star-formation rate of 50 M Q yr" 1 (for a bolometric luminosity 
of 2.8 x 10 11 Lq each) and that Arp 220 is a proton calorime- 
ter, then I find: 

C7 , 4 . 8xl0 -n s ec-.(^)" 2 (^)(^)"'. 

(12) 

Thus, even deep molecular cores in Arp 220's nuclei should 
have an ionization rate comparable to the molecular clouds of 
the Milky Way. 

A further effect which may enhance the GRDR ionization 
rate is dumpiness in the gamma-ray emission of the star- 
burst. I assumed that the gamma-ray emission was evenly dis- 
tributed across the starburst. However, some molecular clouds 
in starbursts may be located near gamma-ray sources. The 
outer layers of the molecular cloud likely themselves gener- 
ate gamma rays as CRs are stopped by them, enhancing the 
interior gamma-ray flux and g amma-ray ionization rate (e.g., 
lUmebavashi & Nakanolll981l) . 

2.2. Comparison with CRDR Ionization Rates 

How does the ionization rate in GRDRs compare with that 
in the surrounding CRDRs? One way of estimating the ioniza- 
tion rate from CRs is to scale the Milky Way molecular cloud 
ionization rate by the energy density of CRs. A potential prob- 
lem with this approach, however, is that the CR spectrum has 
a different shape in starbursts than in the Milky Way. Whereas 
the Milky Way has a E~ 2 1 gamma-ray s pectrum because CR 
protons escape through diffusion (e.g.. iGinzburg & Ptuskinl 
119761) . starbursts are observed to have hard E~ 2 1 gamma-ray 
spectr a from strong advection and/or pionic losses for pro - 
tons (lAcero et al.l2009tlAcciari et al.l2009l:lAbdo et al.l2010h . 
Therefore, some of the larger CR energy density in starbursts 
is in very high energy particles which are not responsible for 
ionizationQ Furthermore, this approach requires a CR ion- 
ization rate for the Milky Way and a Milky Way CR energy 
density. 

Another way to estimate the CR ionization rate is to use the 
injected power in CR protons that is lost to ionization: 

L ion =f^Kf? n (K)dE, (13) 

5 The spectrum of CRs can be approximated by multiplying the injection 
spectrum times a lifetime at each energy, and from this we can get the energy 
density of particles at, say, 1 GeV (EdU jdE = E 2 dN/dE). Then in princi- 
ple the ratio of the starburst and Milky Way EdU / dE at < 1 GeV can be 
used to est imate the ionization ra te. This is effectively done in Such kov et all 
1 1993) and Papadopoulos 1 2010): the CR energy density is related to the star- 
formation rates (injection) divided by an effective escape velocity, which is 
related to the GeV diffusive escape time in the Milky Way and the (fastest ob- 
served) wind speed in starbursts. In any case, however, the calculation does 
require more care than simply computing the integrated energy density. 
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where f(' on (K) is the fraction of energy of a CR with kinetic 
energy K that is lost to ionization and dQ/dE is the injection 
spectrum of CR protons. The fraction of energy going into 
ionization is roughly: 



fi 



P , 

ion 



flife 
* ion 



tr 

1 + J2n. 
U 



-i -1 



^wind 



(14) 



Supposing the pionic losses to have a threshold of K 
140 MeV, the rest energy of a pion, 



fP 

J ion 



-i -1 



1 + _ion_ 

^wind 

t? l 



(K < 140 MeV) 
-i (15) 
(K> 140 MeV). 



The lifetime of protons to ionization losses is approximately 



159 Myrf 



V cm - - 



/ K 



VGeV 







(16) 



where /? = \/ 1 - 1 /7 2 is the speed of the proton divided 
by c ([Torres 2004, and references therein). This gives us, 
roughly, an ionization to pionic lifetime ratio of tf Jt n ~ 
3.18(tf/GeV)/3. 

Now suppose the injection spectrum is a power-law in 
total energy, extending from E = m p c 2 (K = 0) to infinity: 
dQ/dE = Qo(E/m p c 2 )~ p . The normalization of this power 
law is set by L CR = J^ 2 Qo(E -m p c 2 )(E /m p c 2 )~ p dE, giving 
Qo = Lc R [l/(p-2)- l/(p- \)Y\m p c 2 T 2 . When p = 2.2 and 
f„ = 1, I find that L; on = 0.047Lcr. Another commonly as- 
sumed injection spectrum is a power law in momentum. For 
p = 2.2 and f„ = l f with minimum kinetic energies ranging 
from 1 MeV to 100 MeV, I find L ion = (0.10- 0.15)L C r- I 
therefore adopt Lj on ss O.ILcr as a fiducial value in the calcu- 
lations below. 

The CR ionization rate is 



CCR = 



E ion V{n H ) E ion M H 



(17) 



after defining ?7j on = Lj on /LcR. In this equation, Mu is the gas 
mass in the starburst. Scaling to values typical in Arp 220's 
nuclei, I find 



Ccr = 4.4x10" 



' I5 sec" 1 



SFR 



50 M yr 1 / V 4x 10 8 M 



Me 



30 eV 
(18) 



This is ~ 100 times that of Milky Way molecular clouds or an 
Arp 220 GRDR. 

A convenient expression for the ratio of the ionization rates 
in a GRDR to a CRDR is the ratio of equations l2l and [TTl 



C 7 _ f^zfZM n n) _ /^ 7 Z^d( S s) 



Ccr 3?7i n 6rn on m H 

where I have assumed a disk geometry for the volume (V 
2-kR l K). In Arp 220, this ratio reaches 



(19) 



tOMOfJ . n (Ss) _ 2 

Ccr Vl°g cm 



c 7 



f 

•' ion 

o.i 



or 



(20) 



Note this ratio increases even once proton calorimetry is at- 
tained. This is because at high densities, the pionic lifetime 
continues to decrease as n#, so that the energy density of CRs 
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FIG. 3. — Plot of ionization fraction of nu = 10 6 cm" 3 material in GRDRs, 
using the Schmidt Law. The plot also shows the ratio of the ambipolar diffu- 
sion and free-fall times. At high gas surface densities, ambipolar diffusion in 
GRDRs is slower than the free-fall time, suggesting that magnetic fields can 
regulate gas cloud collapse. However, unlike in CRDRs, ?ad / % < 100, so 
the collapse is not slowed far below observations. 

per supernova decreases; the energy density of gamma rays 
per supernova instead remains constant. Only when r 7Z > 1 
for the entire starburst does the gamma-ray flux also start to 
diminish. 

2.3. The Ionization Fraction of GRDRs: Consequences for 
Ambipolar Diffusion 

McKee ( 1989) gives the ionization fraction of a cloud with 
density n H as 



x e = 1 x 10 7 rJ 



»ch 



1/2 



1 + 



»ch 

4n„ 



1/2 



»ch 

4n H 



1/2 



(2D 

Here, n c h ~ 1000r 2 d (C/10 17 sec ') is a characteristic den- 
sity, and r ga - is the g as-to-dust ratio divided by 100 (see also 
Papadopoulos 2010). For most GRDRs, the ionization rate is 
low enough that n cn -C «// and this can be approximated as 



GRDR 



= 3.2x 10" 



1/2 



n u 



10 6 cm" 



-1/2 



(22) 



The ionization fraction when nu = 10 6 cm" 3 is plotted in Fig- 
ure |3] In Arp 220, this fraction rises to x^ RDR « 10" 8 for 
nu = 10 6 cm" 3 . However, in using equation|2T] I have ignored 
the possible appearence of different chemical effects in the 
very low ionization environ ments of GRD Rs, such as charge 
being carried by dust grains (lMcKedll989l) . 

The ionization fraction is important for determining the 
ambipolar diffusion time, the time for the magnetic field to 
slip away from the neutral gas. If the magnetic field re- 
mains strong in a cl oud, free-falling gravitational collapse 
cannot proceed (e.g., Mouschovias & Spitze^ [T976l) . When 
the ambipolar diffusion time is longer than the free-fall time 
and the magnetic fields are strong, the collapse instead pro- 
ceeds quasi-statically o n the ambipolar diffusion ti me (e.g., 
Mestel & Spitzen [l956t see also the discussion in ICrutcher] 
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1999). In most Milky Way starless cores, the ambipolar 
diffusion timescale is roug hly an order of magnitude longer 
than the free-fall time (e.gJCaselli et al.ll99"8tlWilliams et alJ 
119981: iMaret & Berginl 120071: iHezareh et all 120081) . There is 
some recent evidence that at least some of the very deepest 
Milky Way cores may have low ionization fractions and am- 
bipolar dif fusion times approximately equal to the free-fall 
time (e.g. JCaselh et al.ll2002t iBergin & Tafallall2007l) . Given 
that strong magnetic fields are observed b y Zeeman split- 
ting i n dense OH maser regions in ULIRGs (Robishaw et al. 
2008), it seems plausible that magnetic fields regulate star for- 
mation in starbursts, so understanding the ambipolar diffusion 
ra te is important . 

iMcKeel (1989) also gives the ambipolar diffusion time as 
(ad ~ 1-6 x 10 6 (jc,./10~ 8 ) yr. In the low ionization environ- 
ments of GRDRs with n c h -C «h, we have 



,GRDR 



= 5.1 x 10 5 yr 



10" 17 sec" 



1/2 



"// 



10 6 . 



-1/2 



(23) 



We can compare this with the local free-fall time tg 
y/3iv/(32Gm H n H ): 



to = 5.1 x 10 4 yr 



10 6 cnr 



-1/2 



(24) 



The interesting point for star-formation is that both the am- 

bipolar diffusion time and free-fall time scale with n H 1 . Thus 
all GRDRs in a given starburst galaxy will have the same ratio 
of ambipolar and free-fall times, regardless of their overden- 
sity. This ratio is 



fff 



C 7 



1/2 



(25) 



The ratio is also plotted in Figure |3j from the Schmidt Law, 
we expect ambipolar diffusion in GRDRs to be slower than 
free-fall in galaxies with (S g ) > 0.5 g cm" 2 . Even in M82, 
this ratio is ~ 1, so that ambipolar diffusion is no quicker 
than free-fall time. But in a starburst like Arp 220, this ra- 
tio reaches 22. 

The ratio of ambipolar diffusion timescale to free-fall time 
in dense starburst GRDRs is therefore similar to those ob- 
served in Milky Way cores, ~ 1 — 100. This means that 
star-formation may proceed similarly in dense starbursts and 
normal galaxies. By contrast in CRDRs, the predicted am- 
bipolar diffusion time is vastly longer than the free-fall time; 
star-formation would effectively be halted unless some mech- 
anism, perhaps turbulence, speeds up ambipolar diffusion 
(Papadopoulos 2010). It is interesting that in Arp 220, ?ad DR 
approaches ~ 100%, the empirically observed timescale fo r 
star-formation (e.g., Kennicutt 1998; Krum hoTz et al.l 120121) . 
This suggests that in Arp 220, ambipolar diffusion in GRDRs 
may alter the star-formation efficiency slightly. 

2.4. The Temperature of GRDRs 

Papadopo ulos! (120101) gives a convenient expression for the 
minimum gas kinetic temperature of a CRDR which also ap- 
plies for GRDRs, assuming no dust or turbulent heating: 



7; mm = 6.3 K [(0.707n4 /2 C_ 17 + 0.186 2 n|) 1/2 - 



0.1 86^ ] 



S/2-,2/3 



(26) 
The 

temperatures for GRDRs and CRDRs are plotted in Figure 2] 



where n 4 = n H /(\0 4 cm" 3 ) and C-n = C/(10" 17 sec" 1 ). 




FIG. 4. — The minimum gas temperature in GRDRs and CRDRs using the 
Schmidt Law. Black is for nu = 10 4 cm -3 , and grey is for hh = 10 s cirT 3 . 
The solid lines are for GRDRs; dashed are for CRDRs with an energy power- 
law injection spectrum; dotted are for CRDRs with a momentum power-law 
injection spectrum. Individual galaxies may be significantly off the plotted 
relations if they do not follow the Schmidt law. 

Simpler expressions can be found for T™ m in the limit of 
high density and low ionization rate, where a Taylor series 
expansion yields: 



rf m = 9.7Kn 4 



2/3,2/3 
S.-17- 



(27) 



Thus, using the expressions for £ 7 above, I derive approxi- 
mate gas kinetic temperatures: 



T^min 
1 k ' 



0.22 K 



/ n H y 2 /3 / 
V 10 6 cm" 3 / V 



2/3 / (E g ) 



gem 



1.6 



•> ion 

o.i 



2/3 



-2/3 



^30 eV 
«£,)< 0.16 gem" 2 ) 



0.064 K( " h Y 
V10 6 cm" 3 / 



- 2 " ( 



0.93 



f 
J_ ion 

o.i 



F 



2/3 



-2/3 



30 eV 
> 0.16 gem" 2 ) 
(28) 

The minimum temperatures of GRDRs are clearly very low: 
only if (and (Esfr)) is very large and the local density 
is relatively small (< 10 6 cm -3 ) do the temperatures exce ed 1 
Kelvin. For the values of C 7 derived for M82 in section 12.11 
I calculate T™* to be 0.5 K if n H = 10 4 cm" 3 and 0.02 K if 
n H = 10 6 cm -3 . For Arp 220's nuclear starbursts, I find T[" r 
is 9 K if n H 



10 4 cm" 3 and 1.3 K if n H 



10 6 cm" 3 . 



k 

For 



comparison, the CRDR temperature in Arp 220 is ~ 42 K 
for n H = 10 4 cm" 3 (using eqn. [18] for £cr)- It seems that in 
most starbursts GRDRs would be much colder than molecular 
clouds in the Milky Way. 

Of course, the temperatures derived in equation |26]are min- 
imum temperatures derived by assuming dust has a temper- 
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ature of 0, so that the gas is cooled by collisions with dust 
grains. In reality, clouds will be embedded in an intense 
far-infrared radiation field (and the CMB), heating its dust 
grains. If the dust temperature is not cooled significantly by 
dust-grain interactions, if there is no source of dust heating 

- whether light from young stars or absorption of molecular 
line emission - within the GRDR, and if the dust is in thermal 
equilibrium with the external radiation field, there will be no 
net flux in dust thermal emission through the cloud because 
of energy conservation, so that the dust temperature should be 
equal to the external effective temperature. 

For a given dust temperature, dust-grain interactions will 
set a gas kinetic temperature floor where the dust-grain heat- 
ing rate T g -d equals the molecular line cooling rate A^e. 
Papadopo ulosl d2010h quotes 

„ „ , 1n -26/ n H \ 2 {Tk\ 1/2 fT dast -T k \ _ 2 

Yo-d = 2.5 x 10 — — erg cm 

s V10 4 cm" 3 / \KJ \ K J B 

(29) 

and 

*--« x, ^(iP=*r(ira) > -'- f '- fl - 

(30) 

Equating these two rates, I find a temperature floor of 

when the density is low and 7* <C 7d ust , and 

W = r dust -1.7K( I ^)" /2 (^|) 5/2 (32) 

when the density is high and 71 m Tdust- GRDRs with hh > 
40000 cm" 3 will have a kinetic temperature near the dust tem- 
perature when Tdust ~ 40 K. At lower densities, if the heating 
and cooling rates are accurate, the gas will be much colder, 
though given the extreme conditions, the low temperatures 
should be approached with caution. 

Another possible source of heating is the dissipation of tur- 
bulence. Most of the molec ular gas in starbursts exis ts in a 
highly turbulent phase (e.g. jDownes & Solomon|[l998l) . Su- 
personic turbulence rapidly dissipates, pa rtly because the su - 
personic motions generates shocks (e.g.. iStone et al.lfl998l) . 
The cores in molecular clouds and protostars achieve densi- 
ties much higher than most of the mass in the turbulent ISM - 
they are said to be "decoupled" from the surrounding medium 

- although such coll apsed structures can be initially generated 
by turbulence (e.g., Kle ssen et al.l l2000b . It is these clouds 
which have the high colu mns and densities where GRDRs are 
most likely to be present. Bergin & Tafalla (2007) argues that 
in these clouds, turbulence is unimportant, based on the ob- 
served relationship between turbulent l ine width and cloud 
size i n Milky Way molecular clouds ( iMversI 119831: iLarsonl 
120051) . 

The cold temperatures of GRDRs have implications for 
their Jeans mass, which may set the minimum mass of stars 
in starbursts and affect their initial mass function (IMF). The 
Jeans mass is 

M J =(^-) V2 (n Hl n H r^. (33) 
where // = 0.75 x 2 + 0.25 x 4 = 2.5 is the mean particle mass 



in the m olecular ISM, in terms of t he atomic mass of hydro- 
gen inn (Papadopou los et al]|201 lh . In dust-heated gas with 
temperatures near Ta OOI , I find the Jeans mass is 

For a dust temperature of 40 K, the Jeans mass reaches a max- 
imum of 1.8 M Q when n# sa 1.3 x 10 5 cm" 3 . Thus, the Jeans 
mass of GRDRs is similar or perhaps slightly higher than 
that in normal galaxies, unlike in CRD Rs, where it may be 
much higher (^10 M ) in starbursts (c.f. Papadopoulos et al. 
l20TlT) . 

_il conclude that the temperature of GRDRs is more likely set 
°by the dust grains being illuminated by infrared radiation or 
possibly turbulence than gamma-ray heating, except perhaps 
in the lower density regions of the most extreme starbursts. 

3. THE TRANSITION FROM COSMIC RAYS TO GAMMA RAYS: THE 
PERILS OF COSMIC RAY SCATTERING 

At first glance, the ionization rates and equilibrium tem- 
peratures of GRDRs seem surprisingly low compared to pre- 
dictions for CRDRs, given the high gamma-ray luminosities 
of starbursts. In a proton calorimetric starburst, L 7 is about 
one-third of Lqr, so why doesn't equation Q] imply that the 
CRDR heating rate and ionization rate are also low? The 
reason is that galaxies are extreme scattering atmospheres for 
CRs, at least for clouds of sufficiently low CR absorption op- 
tical depth. The flux in equation Q] is not the net flux, which 
is relatively small in galaxies, but the total flux, which is 
much larger. Because of the large scattering optical depth 
of a galaxy, a CR can cross a cloud multiple times: this in- 
creases the odds that it will ionize an atom in it, and therefore 
increases the CR ionization and heating rate. However, the 
disadvantage of this is that a CR has a much greater chance 
of being destroyed by pion production or ionization cooling, 
precisely because it can cross the clouds multiple times. 

The slower CRs traverse a molecular cloud, the quicker that 
the CR flux is attenuated (assuming no CR accelerators within 
the cloud), and the more likely that gamma rays provide the 
ionization. At best, however, CRs can traverse the cloud at 
c, with no deflection by magnetic fields within the cloud. I 
consider the ID version of this case in the Appendix, where 
a molecular cloud with an absorption optical depth to CRs t c 
is embedded in a medium with CR absorptivity a out and scat- 
tering coefficient cr out . Far from the cloud, the CR ionization 
rate is C^. In general, a oal 3> a out : galaxies are scattering 
atmospheres for CRs. I find that even with free travel within 
a molecular cloud, the CR flux is attenuated once the pio- 
nic absorption optical depth of the cloud equals -\/a out /cr out 
(much less than 1). To take a specific example, the radio nu- 
clei of Arp 220 have an average density of 10 4 cm" 3 , so the 
pionic lifetime of CRs within them is ~ 5000 yr according 
to eqn. |4] If the (unknown) scattering mean free path within 
them is Adiff = O.lpc, then a cloud effectively becomes opti- 
cally thick to pion production in this model when 

Sabs ~ v /A diff/(cf 7r )S 7r ~ 0.6 g cm" 2 . (36) 

The rapid attenuation of CRs is shown in Figure [5] The short 
effective attenuation length arises simply because CRs cross 
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FIG. 5. — Plot of CR flux in ID model, as a function of depth z, where 
\z/Z\ < 1 is the location of the overdense cloud. The surrounding ISM is as- 
sumed to be similar to Arp 220. Even clouds with t c ~ 0.01 (~ 0.8 g cm -2 £S 
5 X 10 23 cm -2 ) have much reduced CR fluxes. The CR flux reduction 
needed for gamma ray ionization to dominate is shaded in light grey for 
{E g ) = 10 g citT 2 . Cloud sizes of 0.1 (dotted), 1 (dash-dotted), and 10 pc 
(solid) are shown. 

the cloud many times after being scattered by the surrounding 
galaxy. 

On the other hand, it takes several optical depths before 
the CR ionization rate drops below the gamma-ray ionization 
rate. Equation[l4]is an expression for the optical depth tgrdr 
of a cloud that attenuates the CR flux enough so that gamma- 
ray ionization dominates in the center (F(0)/F(oo) = C7/Ccr)- 
These values are plotted in Figure [6] I find CR absorp- 
tion optical depths of tgrdr ~ 1 are sufficient for gamma- 
ray ionization to become more important than CR ioniza- 
tion. Plugging in values for Arp 220's nuclei from eqn. |20] 
(Ccr/Ct ~ 100), I find gamma-ray ionization starts to domi- 
nate once t c > 0.74. For a typical pionic absorbing column 
of 80 g cm" 2 , this comes out to Sgrdr = 59 g cm" 2 . A cloud 
with this column and density n c will have a mass Mgrdr ~ 
(4/3)vrI]3, RDR /(« 2 m 2 f ) w 1.5 x 10 8 M o (« f /10 6 cm" 3 )" 2 . Tur- 
bulence in the molecular ISM does readily generate transient 
density fluctuations, but simulations sho w that typical colum n 
densities should still be near ~ (Eg) (IOstrikeretalJl200lb . 
which is generally less than 10 g cm" 2 . However, turbu- 
lence can also generate Jeans-unstable, collapsing structures 
such as protostars, and these can easily achieve such columns 
dKlessen et alj [2000). Protostars with n c > 10 10 cm" 3 are 
very likely to be GRDRs in Arp 220's nuclei, since Mgrdr ~ 
1 .5 M Q at these densities, even if they are not shielded by a 
surrounding molecular cloud. 

Moreover, the necessary r c for a GRDR grows only loga- 
rithmically with Ccr/ C-y- When CR scattering in the surround- 
ing medium is slow enough (y/ a out / a ut ^ Ct/Ccr)' tne nec " 
essary optical depth is approximately 



tgrdr ; 



-In 



C7/ CcR 



(37) 
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The larger CR lifetimes of weaker starbursts also means that 
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FIG. 6. — The CR absorption optical depth of a molecular cloud in which 
gamma-ray ionization is more important than CR ionization at its center, us- 
ing the ID model that assumes CRs free stream through the molecular cloud. 
From bottom to top, the lines are for a ut/o" ut = 10~ 6 > 10~ 5 , 10" 4 , 10~ 3 , and 

io- 2 . 

CRs cross a cloud more times before being destroyed outside 
of the cloud: the probability of the CR being destroyed in- 
side of the cloud rather than outside goes up. For M82 and 
NGC 253, with Ccr/C 7 ~ 10000 and a typical CR lifetime 
of ~ 10 5 yr, I find r c > 3.6, or Sgrdr = 290 g cm" 2 , when 
a = (0.1 pc)" 1 . Of course, gamma rays themselves start to 
be attenuated over such columns, but as will be seen in the 
next section, cascade emission from the pair e ± can ensure 
the gamma-ray flux is preserved, to order of magnitude, out 
to these columns. 

Of course, since these calculations are ID, they are only il- 
lustrative. In reality, the ID calculations apply only for clouds 
that are essentially large sheets. In three dimensions, CRs 
have more directions to go that do not intersect clouds. There- 
fore, CRs are more likely to miss small, spherical clouds, and 
would cross it fewer times; the CR flux reduction need not be 
so drastic. 

However, the above calculations are conservative in assum- 
ing that CRs can freely enter into a cloud and cross it at the 
speed of light. If instead diffusion is slow inside a molecular 
cloud, the CR flux within it will be much smaller. The col- 
umn through which CRs can penetrate is then E w Inmn ~ 
( yj Dtiif e /3)nmH. The CRs can only diffuse into clouds with a 
column of less than 



S w 0.038 g cm" 



D 



10 28 . 



1/2 



10 4 cm" 



1/2 



(38) 

Thus, if CRs diffuse slowly in molecular clouds, the GR- 
DRs may encompass large fractions of the molecular material 
within a starburst. Note that hard X-rays are also not atten- 
uated under columns < 2.5 g cm" 2 , so they would contribute 
to ionization at such low column densities. The increase of 
the column with density arises because the effective speed 
at which CRs diffuse increases at smaller scales. Below the 
mean free path, the diffusion approximation breaks down, and 
a treatment similar to the one above, where the CRs are as- 
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sumed to free-stream, must be used. 

The actual propagation of CRs in molecular clouds is 
poorly understood. It is possible t hat CRs essentially free- 
stream through molecular clouds (Cesarskv & Volk 1978; 
Morrill 1982), which could occur if the plasma waves 
that n ormally scatter CRs are quickly damped by colli- 
sions (Kul srud & Cesarskvlll971l) . However, there are rea- 
sons to expect slow diffu sion in or near molecular clouds. 
Skilling & Strong (1976) argued that because molecular 
clouds can absorb CRs, there is a net flux of CRs into the 
cloud; the net flux causes plasma instabilities in the ionized 
material surrounding the cloud, generating waves that scatter 
ionizing CRs away from the cloud. More simply, CRs may 
be excluded from molecular clouds simply because the mag- 
netic field is hi gher, which might lead to a smaller diffusion 
constant (e.g., Gabi ci et all 12007b . Observationally, there is 
radio evidence that diffusio n is slow inside the Galactic Cen- 
ter m olecular cloud Sgr B dProtheroe et al.l l2008 ; Jones et al. 
1201 ll) . The diffusion constant also appears to be small in 
the vicinity of the molecular clouds surrounding the super- 
nova remnants W28 and IC 443, although regions outside the 
molecular clouds are included in th e se measurements (e.g., 
iFujita et al.ll2009b iGabici et al.ll2010b iTorres et al.1l2010h : the 
small diffusion constant could be a result of plasma waves 
generated by the CR flux from the near b y supernova remnants 
in the ionized ISM dFuiita et alj|2010h . ICrocker et ail (1201 lb 
interpret the relative gamma-ray dimness of the Galactic Cen- 
ter as being caused by the exclusion of CRs from molecu- 
lar material, although M82 and NGC 253 seem instead to be 
gamma-ray bright. 

4. THE DEVELOPMENT OF THE PAIR CASCADE: THE LONG 
REACH OF GAMMA RAYS 

Gamma rays do not directly ionize gas; instead they pair 
produce e ± that ionize the gas. The largest uncertainty in the 
estimates of the GRDR ionization rate is the efficiency f^ n 
that the energy in these pair e is converted into ionization of 
the ISM. Pair e ± undergo a cascade, emitting high energy ra- 
diation, some of which may themselves pair produce . At 
low energies, ionization losses (with an energy-independent 
dE/dt) always dominate over the other losses, so low en- 
ergy gamma rays are converted into ionization power rela- 
tively efficiently. Above a few hundred MeV, bremsstrahlung 
losses dominate over ionization losses at all densities, since 
dE/dt oc E for bremsstrahlung. Bremsstrahlung radiation 
from an electron of e nergy E e pro duces gamma rays with 
a typical energy E e /2 (Schlickeiser 2002). In optically thin 
clouds, the gamma rays escape, so that bremsstrahlung losses 
waste the pair energy that could go into ionization. In 
clouds of high enough column, though, the bremsstrahlung 
gamma rays can themselves pair produce e ± and have another 
chance to ionize the ISM. At high energy, synchrotron and 
Inverse Compton (IC) losses with dE/dt oc E 2 are the main 
cooling mechanism for pair e ± . The synchrotron radiation 
from of relevant energies is in the form of radio and in- 
frared emission, which either escapes or is absorbed by dust 
grains. Synchrotron emission therefore does not contribute 
to the ionization of GRDRs. IC radiation is ultraviolet light, 
X-rays, and gamma rays. Ultraviolet and X-rays can directly 
ionize atoms. As with bremsstrahlung, gamma rays either es- 
cape or pair produce depending on the cloud column density. 

Estimating f^ on requires modeling of all of these cooling 
processes. Furthermore, the gamma-ray spectral features of 
starburst galaxies must be taken into account. The gamma 



rays of starbursts are thought to be mostly pionic; because of 
the kinematic threshold of pion production, the pionic spec- 
trum should drop away at energies below 70 MeV in dN/dE 
(e.g. JSteckeril 19701) . or a few hundred MeV in E 2 dN/dE (for 
exam ple, ~ 700- 800 M eV for the Galactic pionic spectrum, 
as in lKamae et alj|2006l) . At high energies, the spectrum will 
be a hard power law. As a result, most of the pair e ± will be 
at high energy, where they will not directly ionize the gas but 
instead cool by radiation. Estimating the ionizing effects of 
the resulting gamma-ray cascade is the subject of this section. 

I make many simplifying assumptions. Once again, I use 
a ID two stream model for the gamma rays. In this case, I 
assume that <? ± neither diffuse nor are advected from the re- 
gions they are created (an on-the-spot approximation), which 
would be the case if D — > in the in the dense environments 
of molecular cloudsQ I can phrase the equation of radiative 
transfer in terms of the column Nh = nnz instead of the phys- 
ical depth z. Each generation i of gamma rays produces pair 
e , which in turn generate the gamma rays of the next gen- 
eration /+ 1. The equation of radiative transfer for these as- 
sumptions is: 



dl' ± 
dN H 



= ± 



-fi 7 z/±- 



QC l + Gi" 1 

^brems 



-'IC 



2n f 



(39) 



where for the primary gamma rays (i = 1), Q\^ ms = Q[q = 
0, and <7 7 z = a 7 z/«# is the total cross section for 7Z pair 
production. The intensity I± is a number intensity per unit 
energy here rather than energy intensity. The 1/2 factor for the 
Q terms takes into account that half of the emitted radiation 
would be emitted towards higher columns (in /+), and half 
towards lower columns (in /_), since the direction of e ± are 
likely to be scrambled by magnetic fields. 

Pair production and the radiation from the pairs is calcu- 
lated using 5-function approximations. For 7Z pair produc- 
tion, the pair e ± from a photon of energy E 1 are assumed to 
all be of energy E e = £ 7 /2: 

Ql = 4a lZ (2E e ) J I i (2E e )dn=4(r + (2E e )+r(2E e ))0L yZ (2E e ). 

(40) 

I then estimate the pair e ± spectrum as 



N' e = Q[, x t m 



(41) 



+ t,M 1 includes losses from 



where t Me = [C + Cms + C<* ^ 'ic 
ionization, bremsstrahlung, synchrotron, and Inverse Comp- 
ton. Defining the loss times for each process as E /(—dE/dt), 
I use the Bethe-Bloch ener gy losses for ionization from 
Strong & Moskalenko (1998), the bremsstrahlung loss es in a 
neutral hydrogen gas from Strong & Moskalenko (1998), and 
the standard Thomso n energy loss formulas for s ynchrotron 
and IC emission from Rvbicki & Lightman (11979b . 

Bremsstrahlung emission is assumed to come out at ener- 
gies E 1 = E e /2: 

4jV](2g 7 ) 



Qi 



brems ' 



^brems 



(42) 



Finally IC radiation is assumed to come out at an energy 
£ 7 = 4/3(E e /m e c 2 ) 2 e dRybicki & Lightmanl 11979b . where 

6 By contrast, calculations of ioniz ation in protostars and protop lanetary 
disks assume that CRs stream at c (e.g., Umebavashi & Nakano 1981). While 
this may be appropriate in these environments, over the much larger scales of 
a molecular cloud, magnetic deflection can be important. 
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eo = 2.7 kT is the energy of photons in the background radi- 
ation field. Since IC spreads the energy of one dex in <? ± 
energy into two dex of upscattered photon energy, I have 
E 2 1 Q[ c =E^N e /(2t lc ): 



(43) 



The volumetric energy input in ionization losses is calcu- 
lated as 

' F N' 

e ^dE e . (44) 

ion 

In addition, I assume that Inverse Compton photons with en- 
ergies between 13.6 eV and 1 MeV are immediately absorbed, 
and contribute all of their energy to ionization. EUV and soft 
X-ray (< 10 keV) photons can directly ionize neutral atoms 
through photoabsorption, though in practice, they can also 
heat the ISM by exciting atoms without ionizing it. Hard X- 
rays lose energy through Compton scattering in large column 
depths (> 10 24 cm" 2 ). The energy lost per scattering is hun- 
dreds of eV to many keV, so these scatterings eject electrons 
from atoms, which then can act as secondary electrons to ion- 
ize other atoms. I calculate the volumetric low energy IC en- 
ergy input as 



le-ic - 



r QicdE^. 



Then the ionization rate is just 



n H E K 



(45) 



(46) 



I applied this method to an input gamma-ray spec trum from 
the fiducial models of Lacki & Thompson (2010) of M82, 
which is detected in gamma rays, and Arp 220's east nu- 
cleus. The western nucleus of Arp 220 wo uld likely pro- 
vide roughly similar results (c.f. Torres 2004J), but it is un- 
clear how much it is powered by an active galactic nucleus 
(e.g., iDownes & Eckartl 120071) . The input gamma-ray spec- 
trum is in terms of unabsorbed volumetric power emitted e; 
I convert it to an intensity by multiplying by a midplane- 
to-edge scale height of h, set to 100 pc for M82 and 50 pc 
for Arp 220's east nucleus: I\(Nh = 0) = eh. I assume the 
modeled column is completely opaque, so that /1(/V™ X ) = 0. 
The spectrum extends down to 1 MeV, and I cut it off at 
10 TeV, since higher energy gam ma rays are likely absorbed 
by 77 pair production processes (Torres 2004] llnouelEOlll 
lLacki & Th ompson 2010). I adopt a density of 10 6 cm -3 and 
a magnetic field strength of 10 milliGauss, similar to those 
found by Zeeman spli tting measurements of dense OH masers 
(Robis haw et al.l 120081) . The radiation field has a blackbody 
spectrum with temperature of 40 K for M82 and 50 K for Arp 
220, but is scaled to a characteristic radiation energy density 
t/rad- M82 is optically thin in the far-infrared, so I calculate 
f/iad = L m /(2nR 2 c), where L TR = 5 9 x 1 10 L is the total 
IR luminosity from San ders et al.l (120031) and R = 250 pc is 
the radius. Arp 220's nuclei are opt i cally thick in the far- 
infrared ("e.g..lDownes & Eckart 2007; Sakamoto et al.ll2008l 
Papadopou los et al. 120101) . so I just use the blackbody energy 
density. 

The resulting ionization rates for a large column of material 
are displayed in Figure [7] I find that, to order of magnitude, 



the gamma-ray ionization rat e at low column depths is sim- 
ilar to that derived in section |3T| using ss 0.1. In M82, 
C 7 ~ 3 x 10~ 19 sec, about three times higher than predicted; 
in Arp 220 East, £ 7 is about 8 x 10~ 17 sec" 1 , about twice the 
estimated value. In reality, these values will vary for GRDRs 
with different densities, magnetic field strengths, and radia- 
tion fields. 

In both M82 and Arp 220 East, the ionization rate is main- 
tained to enormous column densities. In the M82 GRDR, 
the ionization rate remains above 10~ 19 sec -1 to 240 g cm" 2 
(N H = 1 .5 x 10 26 cm -2 ), above 10" 20 sec -1 through 790 g cm -2 
(Nh = 4.7 x 10 26 cm" 2 ), and above the 40 K radioactivity- 
induced ionization rate of 10~ 22 sec" 1 through 2100 g cm" 2 
(Nh = 1.3 x 10 27 cm" 2 ). I find that in the Arp 220 East 
GRDR, gamma rays are capable of sustaining ionization rates 
> 5 x 10~ 17 sec" 1 , the canonical value in Milky Way molecu- 
lar clouds, to columns of 7.0 x 10 25 cm" 2 (1 17 g cm" 2 ). Even 
at columns of 2.9 x 10 26 cm" 2 (490 g cm" 2 ), the ionization 
rate is still ~ 10~ 17 sec" 1 , and it drops to 10" 18 sec" 1 only af- 
ter 6.3 x 10 26 cm" 2 (1050 g cm" 2 ). Thus, even if short-lived 
radioisotopes like 26 Al are well-mixed with the molecular gas, 
sustai ning an io nization rate of 10~ 18 - 10" 17 sec" 1 on their 
own <lLackill2012l> . gamma-rays will be the dominant ioniza- 
tion source to columns of several hundred g cm" 2 in Arp 220 
East. It takes a column of 4400 g cm" 2 (Nh = 2.6 x 10 27 cm" 2 ) 
to reduce £ 7 to less than the 40 K radioactivity induced rate of 
10" 22 sec" 1 in Arp 220 East. The reach of GRDRs is extended 
somewhat by high energy cascade e , which can regenerate 
gamma rays through bremsstrahlung and IC emission (short- 
dashed line). 

I conclude that f^ n « 0.1 is a useful approximation, that 
gamma rays can easily provide significant ionization through 
columns of 100 g cm" 2 , and even with attenuation, they may 
be the dominant source of ionization through columns of sev- 
eral thousand g cm" 2 . 

5. OBSERVATIONAL PROSPECTS 

High energy gamma rays ionize material through the pair 
e ± they produce: the chemistry of gamma-ray ionization will 
therefore be the same as for CR e ionization, and CR ion- 
ization in general. The chemical signatures of cosmic ray 
ionization have been considered recently in several recent pa- 
pers in the conte xt of CRDRs (Meiie rinket all 120061 1201 ll 
iBayet etal]l2011l) . The advent of ALMA, SOFIA, and Her- 
schel allows relevant lines, such as high-/ transitions of CO, 
to be measured. GRDRs probe a different regime of ionization 
and temperature not considered by these papers, and might 
have distinct chemistry. Unlike CRDRs where the dust tem- 
perature is less than the gas temperature, I have argued that 
the gas temperature in GRDRs is less than or equal to the dust 
temperature. The gamma-ray ionization rate of GRDRs may 
make them more similar to typical PDRs, but unlike PDRs, 
the columns in GRDRs are so high that there is essentially no 
UV light whatsoever. 

One feature distinguishing G RDR s from CRDRs is that 
GRDRs should be cold (section El). Whereas CRDRs will 
have gas temperatures of up to 50-150 K, GRDRs will have 
minimum gas temperatures < 10 K even in Arp 220 with- 
out dust heating, and maximum gas temperatures equal to 
the dust temperature. The GRDR gas temperatures could be 
much smaller in less extreme starbursts, < 1 K, depending on 
dust-gas interactions and turbulent heating. The presence of 
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FIG . 7 . — Plot of ionization rate using a prediction for the gamma-ray spectrum of M82 (left) and Arp 220's east nucleus (right), using a ID model and assuming 
Hh = 10 6 cm -3 . Cascade are assumed to remain in place (D = 0). I assume that the power in all Inverse Compton radiation between 13.6 eV and 1 MeV 
instantly ionizes the surrounding gas (dotted line). Total ionization rate (solid), ionization rate from primary gamma-rays (long-dashed line), and ionization rate 
from cascade gamma-rays (short-dashed line) are plotted. 



cold, dense gas would therefore support the idea that there are 
regions CRs cannot penetrate into. 

However, a severe difficulty with actually constraining 
GRDR properties with observations arises because of their 
high column densities. The optical depth of a column of gas 
can be parameterized as 



Mum = Nil I M,/A',,) ( y 



(47) 



Planck observations of local Galactic molecular clouds find 
that r /No w 2.3 x 10" 25 cm 2 at A = 250 fjm in molecular gas, 
assuming j3 = 1.8 (Abergelet al. 2011). This means that the 
optical depth is unity out to 

/ S \ 55 
Ai « 1.1 mm I g _ 2 . (48) 
\ 100 g cm 1 J 

Indeed, the west nucleus of Arp 220 is, if any- 
thing, even worse, being optically thick out to mil- 
limeter wavelengths with (S g ) w 10 g cm" 2 ((Nh) ~ 
10 2S cm" 2 ) (e.g.. iDownes & Eckartl l2007t ISakamoto etai] 
l2008t iPapadopoulos et al]l2010h .lfCRs diffuse very slowly 
in dense starburst molecular gas, then gamma-ray ionization 
may dominate for cloud columns as small as 0.01 g cm" 2 
(section [3}, so that GRDRs can be optically thin for wave- 
lengths longer than ~ 6 fim, though X-ray ionization can 
still be important through these columns. However, if CRs 
free-stream through GRDR clouds, then gamma-ray ioniza- 
tion should only be important at columns ~ 100 g cm" 2 . In 
addition to the dust continuum opacity, molecular lines such 
as CO will have additional opacity. 

If GRDRs are really so optically thick, then, their defin- 
ing feature will be their shadows in molecular lines against 
the brighter, hotter CR-ionized material behind them. Be- 
cause GRDR gas is colder than the CRDR gas, and if GRDRs 



are preferentially in dense cores with low levels of tempera- 
ture, the molecular lines should have different kinematics than 
the CRDR gas, with lower turbulent and thermal line widths, 
making them easier to distinguish. If CRs instead pervade 
the entire starburst, they will heat all of the gas to the same 
temperature, and no molecular line shadows should be visi- 
ble. These molecular line shadows may be detectable in the 
nearest starburst galaxies with ALMA, which will be able to 
make detailed and accurate image s with an angular resolu- 
tion of ~ 0."1 dBrown et al.l l2004). At the distance of M82 
(~ 3.6 Mpc), this translates to 1.7 pc of spatial resolution. 
A cloud of this radius and density of 10 6 cm" 3 has a col- 
umn density of T, g = 9 g cm" 2 (Nh = 5 x 10 24 cm" 2 ; mass of 
5 x 10 5 Mq), large enough to be an interesting test of whether 
CR diffusion is fast or slow in starburst molecular clouds. The 
dust in GRDRs will have the same temperature as that outside 
of them, so there are no similar shadows in dust continuum 
emission. 

A final possibility for directly testing gamma-ray ioniza- 
tion is to look for signs of gamma-ray pair production in 
the radio, by looking for synchrotron emission from the 7Z 
pair e^. There have been studies looking in the Milky Way 
for CR penetration into cores by searching for synch rotron 
frompionic secondary e ± (e.g.. lJones et al.ll2008l.l201 ll) . The 
competing signal is the diffuse radio synchrotron emission 
from the starburst. In starbursts, the diffuse e population is 
likely dominated by secondary e ± produced by the same pio- 
nic loss process that produces most of the gamma rays (e.g., 
Rengarajan 20051 lLacki et al.ll2010l) . It is therefore useful to 
relate both the diffuse synchrotron and synchrotron from GR- 
DRs to the gamma-ray emission of the starburst. For GR- 
DRs with optical depth t 7 z < 1, the volumetric synchrotron 
power is related to the incident gamma-ray flux vF u (GeN) as 
^t ynch (GRDR) « (l/2)^(GeV)r 7Z / s ] nch /(2^), where /7 ynch 
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is the efficiency that gamma-rays are converted into syn- 
chrotron radiation. The factor of 1/2 arises because the 
synchrotron frequency depends on the electron energy (and 
primary gamma-ray energy) squared. Setting ^/v(GeV) 

47rz// 1 ,(GeV) and vfj nch < 
find 



(GRDR) = ^e s / nch (GRDR)^/(4vr), I 



^synch (GRDR) 

v/„(GeV) 



: T 7z/ S ynch/2- 



(49) 



For the diffuse synchrotron and gamma-ray emission from 
starbursts, th e ratio of GHz sy nchrotron to GeV luminosity 
is - 20-40 dLacki et al.ll20llb . If / s ] nch re 1, the ratio can 
be up to ~ 1/2 for GRDRs. Assuming vI u (GeV) within the 
GRDR is equal to its mean value throughout the starburst, the 
potentially larger synchrotron radiation efficiency means that 
the radio intensity towards a GRDR can be up to ~ 10-20 
times larger than on sightlines through the diffuse starburst. 
However, this depends on how pair e ± cool in the GRDR; if 
/synch ^ ^.1, then the radio intensity on the GRDR sightline 
will not be enhanced much. Furthermore, it may be difficult 
in practice to distinguish pair e ± made by gamma rays from 
pionic e ± made by CRs, although synchrotron emission from 
a cloud would still be a sign that some kind of high-energy 
particle is penetrating into it. 

It is therefore difficult to directly test for the presence of 
GRDR ionization, because GRDRs have low ionization rates 
in the first place, they have low heating rates, and are probably 
opaque to infrared lines that could inform us about the chem- 
istry in GRDRs. We may only be able to discover cold regions 
where CRs do not penetrate and note that in gamma-ray de- 
tected starbursts like M82, gamma-ray ionization must exist. 
Fortunately, gamma rays are not deflected by magnetic fields, 
unlike CRs, so there is little question that gamma rays can 
penetrate into molecular clouds, though there remain uncer- 
tainties in f h L, because of the intermediate step of the cooling 
of pair e ± . 

An indirect method of studying GRDRs is through their 
effect on star-formation, particularly the Initial Mass Func- 
tion (IMF). If starburst star-forming regions are CRDRs, the 
gas should be heated to temperatures of ^100 K, rais- 
ing the Jeans mass and forcing the IMF to be top-heavy 
(Papadopoulos et al.1 1201 lh . However, GRDRs are much 
colder, resulting in much smaller Jeans masses, perhaps 
even smaller than in the Milky Way since the gas tempera- 
ture s in GRDRs can be less than 1 K at low densities (sec- 
tion |23|i. As I showed in section l2~4l the GRDR Jeans mass 
is less than 2 M Q for all densities, and is roughly 1 M for 
densities of 10 4 cm -3 and 10 6 cm -3 . Observations of the 
IMF of stellar populations formed in starbursts should con- 
strain whether star-formation occurs primarily in warm CR- 
DRs, cold GRDRs, or some mix of the two. There have 
been conflicting results about what the IMF in starbursts 
is; some studies have suggest ed that starburst e n vironments 
have top-heavy IMFs (e . g.. iRieke et alj [19801; iFiger et alj 
[19991; iSmith & Gallagher! 12001b iMcCradv et all 120031) . but 
others are consistent with a normal IMF in starbursts 
(e.g. IMcCradv et al.1 120031: iKim et alJ 120061: iTacconi et alJ 
2008; see a lso the review by iBastian et al.1 120 lofc On the 
other hand, Ivan Dokkum & Conroyl (120101) recently argued 
that elliptical galaxies have a bottom-heavy IMF (see also 
Ivan Dokkum & Conrovll20Tn) . If confirmed, the evidence for 
a bottom-heavy IMF would support the idea that collapsing 
gas in starbursts passes through a GRDR phase cold enough 



to fragment into small mass protostars. 

6. CONCLUSIONS 

Radiation ionizes media with r re 1 most effectively. At 
low optical depth, the radiation simply escapes; at high optical 
depth, the radiation ionizes the outer layers of the medium, but 
the interior is shielded. In order to guarantee ionization of the 
entire ISM, different kinds of ionizing radiation are needed 
with different penetration depths. In starburst galaxies, even 
CRs may be destroyed by dense molecular gas before they 
can ionize and heat molecular gas at the highest columns. I 
have argued that gamma rays, produced by the interaction of 
CRs with that gas, can provide a guaranteed ionization rate 
through columns of 100 g cm" 2 or more. 

I have calculated the basic properties of the resulting GR- 
DRs: 

• The GRDR ionization rate is small, expected to be of 
order 10" 19 sec" 1 in M82 and 10" 16 sec" 1 in Arp 220's 
radio nuclei (Figure [2). The latter is comparable to the 
Milky Way molecular cloud ionization rate. However, 
gamma-ray ionization still is more important than ion- 
ization from long-lived radioactive nuclei such as 40 K 
in the ISM for most starbursts. 

• The ionization fraction of GRDRs is small, ~ 10" 11 - 
10~ 8 for riH = 10 6 cm" 3 (Figure[3]l. But even these small 
ionization fractions are sufficient to slow down ambipo- 
lar diffusion to timescales comparable to or slower than 
the free-fall tim e of gas in starbursts as dense as M82, if 
the lMcKeel ( 119891) formulas hold in such environments. 
This helps dense gas in starbursts retain magnetic fields, 
suggesting star formation proceeds similarly as in dense 
cores in the Milky Way. On the other hand, the ambipo- 
lar diffusion time is less than 100 times the gas free-fall 
time even in Arp 220, which may solve the problem of 
ambipolar diffusion halting star-formation in CRDRs. 

• GRDRs are very cold, with minimum gas temperatures 
of < 10 K even for Arp 220-like starbursts, and < 1 K in 
starbursts like M82 (Figure 0). Most likely, the heating 
of GRDRs comes from dust-gas interactions or possibly 
turbulence rather than gamma-ray heating, raising their 
gas temperature to the ambient dust temperature. 

• The reason the GRDR ionization rate is so much 
smaller than the CRDR ionization rate (Figure [2]l is 
that CRs are scattered many times in starbursts, while 
gamma rays pass freely through them. The CR flux 
through clouds of low optical depth is therefore in- 
creased by a large factor. However, this poses a dis- 
advantage for CR ionization at large columns: the col- 
umn that CRs must penetrate is much higher than that 
for gamma rays. As a result, CRs may be depleted and 
gamma-ray ionization may take over in clouds of col- 
umn < 100 g cm" 2 in Arp 220 (see Figures [5] and [6] 
for the ID case). If CR diffusion is slow, even clouds 
with column as small as 0.01 g cm" 2 may be GRDRs, 
depending on the presence of X-rays. 

• Calculations of the actual ionization rate £ 7 in GR- 
DRs are complicated by the cascade that develops from 
jZ pair production. I find that assuming = 10% 
as an ionization efficiency works as a rough guide 
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to the ionization rate in densities of 10 6 cm 3 (sec- 
tionHJi. The gamma-ray ionization rate does not appre- 
ciably drop until columns of 100 g cm" 2 . Gamma-ray 
ionization can dominate over 40 K radioactivity (£ 7 > 
10" 22 g cm" 2 ) through columns of several thousand 
g cm" 2 . 

• Observations of GRDRs are likely to be difficult, be- 
cause they occur in regions of columns high enough to 
obscure even the infrared and submillimeter (eqn. |48). 
In principle, the chemical signatures should be similar 
to low levels of cosmic ray ionization. Because they 
are relatively cold, GRDRs would appear as shadows in 
molecular lines against the hotter CRDR regions behind 
them, and would have different line kinematics than 
the background CRDRs. These molecular line shad- 
ows may be visible in nearby starbursts like M82 with 
ALMA. Synchrotron emission from the pair e ± in GR- 
DRs may also be visible in the radio, though it is un- 
clear if such observations could distinguish gamma-ray 
produced from pionic made by CR protons. 

• The cold gas of GRDRs means that their Jeans masses 
can be very small (equation [34i>. If star-forming gas 
passes through a GRDR phase, it can fragment on small 
scales, and may have a normal or even bottom-heavy 
IMF. This is in contrast to CRDRs, which are much hot- 
ter and must have a top-heavy IMF. 

While my focus in this paper has been on gamma rays from 
starbursts, GRDRs can also be present around AGNs. Indeed, 
while gamma ray production in starbursts is limited by the 
relatively small fraction of bolometric luminosity that goes 
into CRs, these limits do not apply to AGNs. Furthermore, 
the i nverse square depen dence of the flux on the distance 
(c.f.. iPapadopou los 2010) actually helps for gas near AGNs. 
Therefore, the gamma-ray ionization rate near AGNs may be 
much higher, and GRDRs may set the conditions for star for- 
mation that appears to happen near gam ma-ray bright AGNs 
(Paumard et alJl2006t iDavies et al.ll2007h . 

GRDRs are the next rung in the hierarchy of ionization re- 
gions after CRDRs (Figure [T). It is natural to ask whether 
anything is even more penetrating than gamma rays. One 
known source of ionization at arbitrarily high columns is ra- 
dioactivity from unstable isotopes in the ISM, such as 40 K 
(lCameron|[T96l lUmebavashi & Nakanol fl98lt) . The 26 Al 
produced by the many young stars in starbursts may in- 
crease the radioactive ioniza tion rate of the ISM further 
(lUmebayashi & Nakanoll2009i) . If 26 Al is rapidly mixed with 
the molecular gas of starbursts, it should sustain ionization 
rates of 10" I8 -10" 17 sec" 1 dLackill2012l) . Then 26 Al would 
dominate the ionization of dense molecular gas in weaker star- 
bursts (including M82) and gamma rays would dominate the 
ionization rate in stronger starbursts (like Arp 220's nuclei), 



as shown in Figure [2] Even if 26 Al does not contribute to 
the ionization rate, the gamma rays provide a relatively cer- 
tain source of ionization in extremely dense gas. As noted in 
section|U gamma rays can provide greater ionizing rates than 
40 K radioactivity even to column depths of several thousand 
g cm" 2 in starbursts. 

This would be very relevant for protostars and proto- 
planetary disks, which can achieve column densities of 
1000 g cm" 2 . It is thought that protoplanetary disks in the 
Milky Way have "dead zones" where there is essentially no 
ionization at column densities greater than 100 g cm" 2 . In the 
dead zone, the gas is no longer affected by ma gnetic fields 
an d accretion ma y be halted (e.g., Gammie 1996). According 
to iGammid (119961) . dead zones appear in regions with x e < 
10" 13 ; for regions of characteristic density of = 10 13 cm" 3 , 
this corresponds to an ionization rate of £ rj 5 x 10~ 20 sec" 1 , 
smaller than the gamma-ray ionization rate in Arp 220 by sev- 
eral orders of magnitude. Hence, in starbursts the dead zones 
may be reduced or even non-existent because of gamma-ray 
ionizatio n, which may affect the process of planet forma- 
tion (c.f. iFatuzzo et al J 12006. which discusses the effects of 
high CR ionization on protoplanetary disks). Another possi- 
ble region where gamma rays would be attenuated are nuclear 
torii surrounding AGNs, which can have clumps with column 
densities approachin g Nh = 10 27 cm" 2 (S g = 1700 g cm" 2 ) 
(Hopk ins et al.ll2012l and references therein). 

A more speculative radiation source is dark matter annihi- 
lation or decay. Dark matter is distributed throughout galax- 
ies and is not stopped by matter of any reasonable density, 
although it would be concentrated more towards the galactic 
center, where any active nucleus would be. High concentra- 
tions of dark matter may produce high energy radiation which 
will cascade down and contribute to ionization. The most 
extreme regions of starbursts or active nuclei would then be 
Dark Matter Dominated Regions (DMDRs), and indeed the 
entire early Universe before reionization may have been a gi- 
ant DMDR, when there was no other source of ionizing ra- 
diatio n (e.g..|Padmanabhan & Fink beinerl2005l: iMapelli et all 
2006). 
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APPENDIX 

I consider ID radiative transfer in a proton calorimetric starburst. For the ID problem, there are just two sightlines, one 
extending to +z with an intensity 1+ and the other extending to — z with an intensity /_. The equation of radiative transfer is 

d ^ =± ^m- a(z)I± - a(z)(I± -j^ (i) 



where J = (7+ + /_)/2 is the angle-averaged intensity at that z, a(z) is the inverse of the mean free path for CR destruction 
(absorption coefficient), and a(z) is the inverse of the mean free path for CR scattering, and Q(z) is the rate that CRs are injected 
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into the ISM dRvbicki & Lightman ll979t) . Suppose a molecular cloud sits at -Z < z < Z, through which CRs free stream at c. 
Outside of the cloud, CRs can be either scattered or absorbed; I assume this surrounding, large-scale medium is homogeneous. 
The surrounding starburst is also assumed here to be infinitely large, which should yield the correct behavior as long as the 
starburst is large compared to the typical diffusion length of protons. The cloud, being denser than the surrounding medium has 
a different absorption coefficient than the surrounding medium. Furthermore, I assume that no CRs are emitted inside the cloud. 
Therefore: 



Q(z)' 



a c (z<Z) 
a out (z > Z) 

(z < Z) 

Cr ou t (Z > Z) 

(z < Z) 

gout (Z>Z) 



(2) 
(3) 
(4) 



a r Z. 



I define the CR absorption optical depth from the edge of the cloud to its center (z = Z to z = 0) as r c : 
The equation can be solved for z > Z by differentiating both sides by z: 

-jY = -\®om+ (Tout/ 2)— + <r ut/ 2 — 

and then using Q]to find, in terms of I+, substitutions for dl-/dz and then /_ = 2 / a ont [dI + / dz - Q ou t/2 + ((Xout+ &out/2)I+]- The final 
differential equation is then 
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After requiring I± to be finite as z — > oo, the solution to the differential equation for z > Z is: 
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(7) 
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where I define /i = ^\ + cr out /a out . In these equations, C is a constant set by the boundary conditions at the edge of the molecular 
cloud (z = Z). From symmetry, the intensity of CRs entering into the cloud at z = +Z must equal the intensity of CRs entering 
into the cloud at z = -Z, so that I-(Z) = I+(-Z). In addition, since there is no CR scattering within the cloud, the CR intensity 
emerging out of the cloud is just the intensity emerging into the cloud attenuated by absorption: /+(Z) = I-(Z)e~ 2 ° cZ and likewise 
7_(-Z) = I + (-Z)e~ 2a ' z . To find C, I therefore impose the additional condition that /+(Z) = 7_(Z)e~ 2Tt . This gives me for z > Z 
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(9) 
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Now I can compare the CR flux F = J Idil = /+ + /_ inside the cloud to the CR flux far away from the cloud. At the center of 
the cloud, z = 0, the CR intensity is the intensity entering into the cloud after being attenuated: /+(0) = /-(0) = I+(Z)e Tc = 7_(Z)e" Tl . 
Far from the cloud, as z — > oo, /+(oo) = /-(oo) = <2out/(2awX which simply means that the CR density is equal to their injection 
rate times their lifetime. Taking the ratio of the fluxes, I have 



F(Q) 
F(oo) 



= exp(-r c ) 



1 



cr out / (2aout) - ( 1 + CT ut/(2a ut) - A*)e 2t 



(11) 



In a galaxy, the CR scattering mean free path is thought to be much smaller than the absorbing mean free path, so that a out 3> a out . 
This ratio is roughly 

m) - ' (12) 



F(oo) V CT °ut 1 - e" 2 ^ + 2 ^a om /a out e- 



As it turns out, even clouds that are optically thin to CR absorption on one pass can suppress the CR flux within them. If 
t c <c 1 , 1 have 

F(0) - (13) 



F(oo) Y CT out ^ow /cr ut + ( 1 - 2 ^a out / C ut )t c 

For very small cloud absorption optical depths, the flux within the cloud is the same at infinity. However, once r c reaches 
~ \J Q^out / Tout <C 1, the flux begins to drop inversely with r c . This is despite the fact that CRs free stream within the cloud in this 
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model and are not trapped inside it: the flux is low because they can make multiple passes through a cloud and have multiple 
chances to be destroyed. 

Equation Q~2] lets me solve for the approximate column (in the 1-D planar case) when gamma-ray ionization dominates over 
cosmic ray ionization. Setting F(0)/F(oo) = C7/CCR where Ccr is the CR ionization rate outside of the cloud, I get 



tgrdr : 



-In 



(Ccr/ C7)\/aout7^out+ y aout/fTout(CcR/ C7) 2 + l-2 y / a oul / 



fout 



l-2^a out /c 



(14) 



This calculation neglects absorption of gamma rays within the cloud, which will start to matter over columns of ~ 200 g cm 2 
(tgrdr ^ !)• In the limit that scattering outside the cloud is slow and the exterior CR ionization rate is much greater than the 
gamma-ray ionization rate, with y/a 0Ut /(T 0Ut <C C7/CCR' equation [14] can be simplified to 

(C7ASD 



tgrdr : 



-In 



2-\/«out/o'out 



(15) 
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